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We report the discovery of a Compton-thick AGN and of intense star-formation activity in the nucleus and disk, respectively, of 
the nearly edge-on superwind galaxy NGC 4666. Spatially unresolved emission is detected by BeppoSAX only at energies < 10 
keV, whereas spatially resolved emission from the whole disk is detected by XMM-Newton. A prominent (EW ~ 1 — 2 keV) 
vQ . emission line at ~6.4 keV is detected by both instruments. From the XMM-Newton data alone the line is spectrally localized 

f^ ' at i? ~ 6.42 ± 0.03 keV, and seems to be spatially concentrated in the nuclear region of NGC 4666. This, together with the 

^^ , presence of a flat (F ~ 1.3) continuum in the nuclear region, suggests the existence of a strongly absorbed (i.e., Compton-thick) 

^^ ■ AGN, whose intrinsic 2-10 keV luminosity is estimated to be I/2-io^ 2 x 10*^ erg s^^. At energies ^ I keV the integrated 

{BeppoSAX) spectrum is dominated by a ~0.25 keV thermal gas component distributed throughout the disk (resolved by XMM- 
Newton). At energies ~2-I0 keV, the integrated spectrum is dominated by a steep (F^ 2) power-law (PL) component. The 
latter emission is likely due to unresolved sources with luminosity L ~ 10"^* — 10^^ erg s~^ that are most likely accreting 
O I binaries (with BH masses <?>Mq). Such binaries, which are known to dominate the X-ray point-source luminosity in nearby 

star-forming galaxies, have F ~ 2 PL spectra in the relevant energy range. A F ~ 1.8 PL contribution from Compton scattering 
^ ' of (the radio-emitting) relativistic electrons by the ambient FIR photons may add a truly diffuse component to the 2-10 keV 

emission. 
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1 . Introduction (a) The torus, which is presumed to be essentialy neutral, gives 

rise to the 'Compton reflection continuum' or 'cold reflection' 

The Compton-thick subclass of Seyfert 2 (Sy2) galaxies con- ^^^^ q^^^^^ ^ p^^i^^ I991) that can be explained as follows. 

sists of objects for which the column density of absorbing mat- ^^^^^^^ ^„ -^^^^^ power-law (PL) continuum iiTadiates some 

ter along the l.o.s. to the circumnuclear torus is Nn > ctt - 'cold' (i.e., T< lO^ °K) gas with solar chemical abundance. 

1.5 X 10^4 ^^ 1 A column of a few times this critical value ^^^ ^^^^cis^ X-ray spectrum contains the Fe-Ka (6.40 keV), 

IS sufficient to strongly depress the transmitted mtensity also m p^.^^ (7 q^ keV), Ni-Ka (7.48 keV) fluorescence lines, along 

hard X-rays, since after a few scatterings the radiation is red- ^-^^^ ^^e Fe-K-shell absorption edge at ^7. 11 keV. Due to the 

shifted into the photoelectrically-dominated regime. ^^^^ abundance and fluorescence yield of iron, a significant 

However, even for very high values of A^h, nuclear radia- fraction of the incident radiation at energies E > 7.11 keV 

tion can be observed by the radiation scattered from ffl) the vis- yields Fe-Ka fluorescent line at ^ = 6.40 keV. Compton 

ible part of the inner surface of the circum-nuclear torus (that down-scattering within the gas causes a fall-off in the contin- 

blocks direct X-ray emission entirely from our view; Antonucci uum at high energies, a smearing of the Fe-K-shell absorption 

1993), and from scattering by (b) the warm medium that scat- edge, and the formation of low-energy wings of the fluores- 

ters and polarizes the optical broad lines (Antonucci & Miller cence lines. The competition within the slab between multiple 

1985). The two scattered components are spectrally different. electron down-scattering of high energy photons and photo- 

electric absorption of low-energy photons leads to a 'Compton 
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reflection hump' in the reflected X-ray continuum in the ^ 20- 
100 keV band. 

(b) The warm mirror, which is presumed to be highly ionized 
(i.e., T^ 10^ °K) and of low optical thickness, is located above 
the obscuring torus and electron-scatters the continuum pho- 
tons coming from the central source (Krolik & Kallman 1987). 
If the ionized scattering medium is optically thin, photoelectric 
absorption is of little importance on the electron-scattered con- 
tinuum, so the emergent spectrum essentially retains the spec- 
tral shape of the primary emission (except for a high-energy 
cutoff due to Compton down-scattering). This same warm scat- 
tering medium most likely produces the higher-energy Fe-Ka 
lines from, respectively. He-like iron (Fe XXV) at ^6.7 and 
H-hke iron (Fe XXVI) at 6.97 keV. 

A combination of warm and cold reflections is, there- 
fore, a plausible model for the observed X-ray spectra of Sy2 
sources hidden in galactic nuclei (Matt et al. 1996, 1997). The 
relative strength of the cold and warm reflections can pos- 
sibly be determined from fitting the composite model to 2- 
10 keV ASCA measurements of the proto-typical Sy2 galaxy 
NGC 1068. From the fraction of 2-10 keV flux that comes 
from cold reflection, and an independent estimate of the in- 
trinsic L^_io - 7 X 10^3 erg s^^ from the [O I1I]A5007 lu- 
minosity, and a 2-10 keV albedo of 0.022 for the whole inner 
wall of the torus, Iwasawa et al. (1997) deduced that a fraction 
/c — 4.4 X 10~^(L43 g)^^ of the inner surface of the torus was 
visible to us. Such a small fraction of the reflecting inner sur- 
face of the torus indicates that the torus is viewed nearly edge- 
on. Iwasawa et al. (1997) further deduced the warm-scattered 
luminosity fraction to be /w — 7.8 x 10~''(L43 g)~^. The 
warm- to-cold reflection luminosity ratio in NGC 1068 is then 
0.018. For a less extreme inclination of the absorbing torus - 
i.e. for a higher fraction of the reflecting surface -the 2-10 keV 
spectrum of a Sy2-like source is even more dominated by cold 
reflection. Thus, we may conclude that the spectral features of 
cold reflection, i.e. a 'flat' (F ~ 1.3) 2-10 keV continuum and 
a 6.4 keV Ka line from 'cold' iron, constitute the X-ray sig- 
nature of a strongly absorbed (i.e., Compton-thick) Sy2-like 
source. Based on these results, a strategy aimed at unveiling 
strongly absorbed Active Galactic Nuclei (AGNs) harbored in 
galactic nuclei is to make spectral observations of galaxies in 
the band 2-10 keV, and to look for a prominent 6.4 keV line 
superimposed on a rather flat underlying continuum. 

A substantial link between circum-nuclear starbursts (SBs) 
and AGNs may be established if both occur in the same galaxy. 
If the SB-driven turbulence in the gas increases the mass in- 
flow rate onto a central black hole, a higher emission from the 
central engine will ensue (e.g., Veilleux 2001). Since a circum- 
nuclear SB is located farther out than the torus and its obscur- 
ing material is distributed quasi-isotropically as seen from the 
central source, the nuclear source will be further obscured by 
the SB so that the emerging (cold-reflected) X-ray spectrum 
will be further hardened by SB-related photoelectric absorp- 
tion (Fabian et al. 1998). Roughly 50% of all Sy2 galaxies con- 
tain circumnuclear SBs (Gonzales Delgado et al. 2001): these 
composite objects do tend to be more obscured than those lack- 
ing SBs. It is indeed reasonable to conjecture that such extra 



obscuration is caused by the SB, because star-formation rates 
(SFRs) per unit area of --50-100 Mq yr^^ kpc"^ (typical of 
the innermost 0. 1 kpc in such composite galaxies; see Gonzales 
Delgado et al. 1998) imply A^h ~ 10^^ cm^^ (see Kennicutt 
1998). A circum-nuclear SB can also be physically connected 
to a torus. If the mass inflow rate through the accretion disk is 
higher than the rate of power dissipation in the disk, its opac- 
ity will increase, and the disk will puff up into a geometrically 
thick torus with a large covering angle (as seen from the central 
engine). This will quickly block the primary AGN continuum 
and will result in emission of the fluorescent Fe-Ka line and 
flat reflection continuum. Indeed, among the heavily obscured 
(A^H > 10^"' cm^^) sources studied by Levenson et al. (2002), 
those that exhibit the largest equivalent widths (EWs) have con- 
centrated circumnuclear SBs, whereas in the only galaxy cer- 
tainly lacking a SB, Mrk 3, the Fe-Ka line is relatively weak 
(EW - 0.7). 

Well-known examples of coexisting spectral signatures of 
star formation and Compton-thick AGNs include cases of both 
type-1 and type-2 view of the central AGN. (These views de- 
fine Syl and Sy2 morphology in the 'Unified AGN Model', 
e.g. Urry & Padovani 1995.) Examples of the former type in- 
clude Arp 299, NGC 6240, and NGC 4945: in these objects 
the absorbing matter is moderately thick to Compton scattering 
(A^H ^ a few X 10^^ cm^^), so the the primary AGN spectrum 
emerges at (source restframe) energies E^ 10 keV, within the 
BeppoSAX observing window (Delia Ceca et al. 2002 and Ballo 
et al. 2004; Vignati et al. 1999; Guainazzi et al. 2000). The 
archetypal example of the latter type is NGC 1068 for which 
BeppoSAX data support a model envisaging a mixture of both 
cold and warm reflections of an otherwise unseen primary nu- 
clear continuum, the cold reflection component being the dom- 
inant one in the 20-100 keV band (Matt et al. 1997). 

NGC 4666 is a starburst galaxy (SBG) seen nearly edge-on 
(see Table 1). Based on multi-frequency observations, Dahlem 
et al. (1997) have detected an outflow cone - associated with 
a galactic superwind - emanating from a central SB of ~3.3 
kpc in radius, and having an opening angle of 30° ± 10°. 
The outflow is traced up to ^7.5 kpc above the disk plane by 
optical emission line filaments, nonthermal radio continuum 
emissiom, and soft X-ray emission from hot gas. The multi- 
wavelength evidence for enhanced star-formation activity in 
NGC 4666 can be summarized as follows: 

(/) Optical spectroscopy exhibits the kinematic signature of an 
outflow, and the observed line ratios indicate shock-heating 
[e.g., following SB-related supernova (SN) explosions] as the 
most likely excitation mechanism for the emission-line gas. 

(//) The existence of a synchrotron radio halo (Sukumar et 
al. 1988; Dahlem et al. 1997) is strong evidence of ongoing 
or very recent SB activity. The radio (energy) index ranges 
from ar ^ 0.8, measured at small galactocentric radii, to 
a^ 1.3, measured at large radii. This is proof of the pres- 
ence, close to the disk, of ongoing injection of fresh rela- 
tivistic electrons whose lifetime against radiative losses, t ~ 
lC^^^{E/GcY)-^{B/^lG)-^ yr, is i - 5 - 50 Myr [for B - 
15 mG (Dahlem et al. 1997) and S - 1 - 10 GeV], i.e. com- 
parable or shorter than a typical SB age (isB ^ 10^ yr). The 



degree of steepening is just what is expected when the electrons 
lose energy radiatively. 

(Hi) At energies <2.4 keV NGC 4666 was previously observed 
with ROSAT (Dahlem et al. 1997, 1998); the 4"-resolution HRI 
image showed faint diffuse emission extending along the ma- 
jor axis and confined mostly to the disk plane, and two point 
sources located ^ 1' along the minor axes above and below 
the nuclear region. The spatially integrated PSPC spectrum was 
fit with a two-component model (Dahlem et al. 1998) consist- 
ing of a PL H- thermal model, with the PL photon index set at 
r = L9, and a derived plasma temperature (with the chemical 
abundances set at the solar value) of kT ~ 0.3llo'52 l^^^- 

In this paper we present the first broad-band (0.1-50 
keV) spectral observations of NGC 4666, obtained using the 
BeppoSAX and XMM-Newton orbiting telescopes. In Sections 
2 and 3 we describe the data reduction and analysis; in Section 
4 we describe the spectral fitting, whose implications are dis- 
cussed in Sections 5 and 6. Our main results are summarized 
in Section 7. Unless explicitly quoted, all the errors reported in 
this paper are at the 90% confidence level for one interesting 
parameter (Ax^ ^ 2.71). 

2. Observations and Data reduction 

2.1. BeppoSAX measurements 

NGC 4666 was observed by BeppoSAX with the three narrow- 
field instruments (NFIs) (see Tables 2 and 3): the Low Energy 
Concentrator Spectrometer (LLCS: Parmar et al. 1997), the 
Medium Energy Concentrator Spectrometer (MECS: Boella et 
al. 1997), and the Phoswich Detector System (PDS: Frontera et 
al. 1997). 

The cleaned and linearized data produced by the BeppoSAX 
Scientific Data Center ' (ASDC) were analyzed using stan- 
dard software packages (XSELECT vl.4, FTOOLS v4.2, and 
XSPEC vl 1.2). For the MECS, we used the event file made by 
merging the data of the two, properly equalized, MECS units. 
At the spatial resolution of the NFI instruments, NGC 4666 was 
not resolved. No significant source flux variation was detected 
over the observing period. 

2.2. XMM-Newton measurements 

NGC 4666 was observed by XMM-Newton on 2002 June 27- 
28, for a total good exposure of ~50 ks with the European 
Photon Imaging Camera (EPIC) detectors. Data were taken 
with the EPIC PN detector (Struder et al. 2001) in extended 
full-frame mode, and with the EPIC metal oxide semiconductor 
(MOS) detectors (Turner et al. 2001) in full-window mode (see 
Tables 2 and 3). The data were reduced using version 5.4.1 of 
the XMM-SAS software, using the standard processing scripts 
(emproc and epproc). 10-12 keV light curves were inspected 
and found to be at a nominal rate for the whole observation. 
Data were selected using event patterns 0-12 (for the MOS) and 
pattern 0-4 (for the PN), and only good X-ray events (using the 
selection expression "FLAG=0" in evselect) were included. 



See_http://www.asdc.asi.it/bepposax| 



Images obtained with the PN detector and with the MOS 1 
are shown in Fig.l. These clearly show that NGC 4666 is ex- 
tended in both energy bands. The extension is most promi- 
nent at the lower energy band, and seem to originate in a truly 
diffuse emission. At higher energies, the emission is primar- 
ily from point-sources, including a nuclear source that is the 
brightest. 



3. Spectral analysis 

3.1. BeppoSAX data 

To maximize the statistics and the s/n ratio, the LECS and 
MECS source counts were extracted from a circular region of 4' 
radius. Background counts were extracted from high Galactic 
latitude "blank" fields (provided by the BeppoSAX ASDC) us- 
ing an extraction region corresponding in size and detector po- 
sition to that used for the source. The PDS spectra extracted 
with the standard pipeline (with the rise-time correction ap- 
plied) provided directly by the BeppoSAX ASDC; the simulta- 
neously measured off-source background was used. Significant 
emission of NGC 4666 was detected at significant levels by the 
LECS and MECS, but not by the PDS instrument (see Table 3). 
Spectral channels corresponding to energies 0.1-4 keV and 
1.8-10 keV have been used respectively for the analysis of 
the LECS and MECS data [as suggested by the BeppoSAX 
'Cookbook' (see Fiore et al. 1999)]; PDS data are in the 13-60 
keV spectral channels. LECS and MECS source counts were 
rebinned to have s/n > 3 in each energy bin; PDS source 
counts were rebinned to have s/n > 2. Standard calibration 
files were used in the fitting procedure; the LECS/MECS and 
PDS/MECS cross-constants normalizations were allowed to 
vary in the ranges proposed by the BeppoSAX 'Cookbook' . 

3.2. XMM-Newton data 

The nuclear spectrum was extracted (for both MOS and PN 
detectors) from a circular region of 15" radius around the nu- 
cleus (shown in Fig.l). Background spectra were extracted 
from offset circles, close to NGC 4666 but free of any back- 
ground source. Standard response matrices and ancillary files 
were used. The background-subtracted spectrum was fitted us- 
ing XSPEC version vll.2, including 0.3-10 keV data for the 
MOS and 0.5-10 keV data for the PN. Data were binned so 
as to have at least 20 counts/bin to enable the use of x^ mini- 
mization during the spectral fits. The plot in Fig. 2 was obtained 
using s/n < 3 (per bin). 



4. I\1ain spectral components 

The basic model components used in various combinations to 
fit BeppoSAX data are a PL, a gaussian line, and the Mekal 
model describing thermal emission from optically thin hot 
plasma (see the XSPEC package). We include photoelectric ab- 
sorption due to Galactic foreground, corresponding to HI col- 
umn density A'n.Gai ~ 1-7 x lO^*' cm^^. Also included is 
instrinsic absorption in NGC 4666 to be determined from the 
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Table 1. General characteristics of NGC 4666. 



Type*"' 


RA('') 
(J2000) 


deg'*") 

(J2000) 


(Mpc) 


^i,0(d) 

(mag) 


logL«(^' 


logiV^-''' 
(cm-2) 




inclC^) 
(deg) 


„i,0(i) 
■f^25 

(kpc) 


logia5-2 
(erg/s) 


log^^'Jio 
(erg/s) 


SABc 


12''45'"0r 


-00°27'38" 


26.3 


10.68 


10.76 


20.23 


10.69 


80° 


17.5 


40.20 


40.40 



^'^' Morphological type, from RC3 (de Vaucouleurs et al. 1991). 

<*^ RA(2000) and DEG(2000), from RC3 (de Vaucouleurs et al. 1991). 

*"' Distance, from Dahlem et al. (1997). 

''*' Total B-magnitude, corrected for Galactic and internal extinction, from RC3 (de Vaucouleurs et al. 1991). 

'■^^ B-band luminosity, derived using the quoted distance and B-magnitude. 

" ^ Galactic HI column density (Dickey & Lockman 1990). 

'*' Far-infrared (FIR) luminosity (from Dahlem et al. 1997). The FIR flux is defined (Helou et al. 1985) as a combination of the IRAS GO/im 

and lOOfim fluxes according to /fir = 1.26 x 10^^^ (2.58 /eo + /loo) erg s^^ cm^^, where /go and /loo are expressed in Jy. 

^^' Inclination angle with respect to line of sight (Dahlem et al. 1998). 

*■*' 25--B-mag photometric radius, corrected for inclination and Galactic absorption, derived using the angular value given by RC3 (de 

Vaucouleurs et al. 1991) and the quoted distance. 

"' 0.5-2 keV luminosity, derived using the quoted distance and the data presented in this paper. Not corrected for Galactic absorption. 

'*' 2-10 keV luminosity, derived using the quoted distance and the data presented in this paper. Not corrected for Galactic absorption. 

Table 2. BeppoSAX and XMM-Newton instrumental parameters. 





FOV" 






Ei-Ei' 




AE^ 






Sensitivity 


d 




SAX 


LEGS 

37' 


MEGS 

56' 


PDS 

~80' 


LEGS MEGS 
0.1^.5 1.6-10 


PDS 

13-300 


LEGS 

0.24 


MEGS 

0.24 


PDS 

9 


LEGS 

10-^^ 


MEGS 

10-" 


PDS 

2x10-'" 


XMM 


MOSl/2 
30' 


PN 

30' 




MOSl/2 PN 
0.3-10 0.5-10 




MOSl/2 
0.13 


PN 
0.13 




MOSl/2 

7xl0-"-' 


PN 

7xl0-''' 


- 



"^ Diameter of the field-of-view. 

b 



Instrument energy range, in keV. 
" Energy resolution, in keV (at 6 keV for LEGS, MEGS, MOSl/2, PN; at 60 keV for PDS). 
"^ Flux from point source after 100 ks exposure (erg cm^^ s^^). 

Table 3. BeppoSAX and XMM-Newton observation log. 





Date 


Exposure 


time " (ks) 


Gount Rats 


;''(10-^cts-^) 


SAX 


2001 Dec 14-18 


LEGS 

42.92 


MEGS PDS 

131.83 58.11 


LEGS 

2.23 


MEGS PDS 

3.76 < 18.0 


XMM 


2002 June 27-28 


MOSl/2 

57.1 


PN 

48.3 


MOSl/2 
~ 0.013 


PN 

0.04 



" On-source net exposure time. The LEGS exposure time is considerably shorter than that of the MEGS, because the LEGS can operate only 

when the spacecraft is not illuminated by the Sun. 

*" Background-subtracted source count rates, with photon counting statistics errors. For the PDS the 90% upper limit is given. 

spectral analysis. The results of the spectral analysis are sum- keV suggest the presence of a soft thermal component, and the 
marized in Table 4. feature at ^^6.5 keV suggests Fe-Kct line emission. 

Multi-component models including different combinations 

Single-component models do not provide an adequate de- of thermal, PL, and gaussian-line components can then be fit- 

scription of the BeppoSAX broad-band spectrum of NGC 4666. ted to the data. A PL-ngaussian-line model (model C), although 

A single unabsorbed PL model (model A) and a single- clearly marking an improvement over models A and B, is still 

temperature thermal model (model B) are rejected at >99.9% unsatisfactory because it does not remove the excess residu- 

confidence level. In Fig.2-left we show the ratio between the als at ^0.5 keV (Fig.2-right). A simple model that provides an 

data and model A. A broad excess is evident at ^0.5 keV, and acceptable statistical fit (model D) includes: i) an unabsorbed 

a line-like excess is seen at ~6.5 keV. The residuals around 0.5 thermal plasma with kT ~ 0.23 ± 0.05 keV (with the chemical 



Fig. 2. Left: residuals of the fit of BeppoS AX data to model A 
(PL with r == 2.15). Right: residuals of the fit of BeppoSAX 
data to model C (PL with T = 2.27, plus a gaussian line at 
E = 6.46 keV): the excess at ^6.5 keV has disappeared, but 
not that at < 1 keV. 



Fig. 1. Soft {top) and hard (bottom) X-ray contours obtained 
with the MOS 1 overplotted on the optical image from the DSS. 
Contours have been obtained after gaussian smoothing. They 
span a surface brightness range of ^ 0.38-9.4 c arcsec^^ in the 
soft band, and ^ 0.3-3.6 c arcsec^^ in the hard band. In the 
lower panel, the 15"-radius annulus marks the nuclear extrac- 
tion region, data from which are shown in Fig.5. 

Table 4. Spectral analysis of BeppoSAX data 



Model 



Parameters 



A: pi 

r 



xVi/ = 49.62/43 
2.15 ±0.15 



B: mekal 
kT (keV) 

Z{Zq) 



X 



Vi/ = 50.24/43 



4.09 ±0.77 
1.0(f) 



C: pl+gauss 

r 

Line: (je; E(keV); EW (keV) 



xVi^ = 38.63/41 

2.27 ±0.17 

0.0(f); 6.46 ±0.55; 1.87 ±0.60 



D: mekal+pl+gauss 
kT (keV) 
Z{Zq) 

r 

Line: a^\ E(keV); EW (keV) 



xVi^ = 33.46/39 

0.23 ±0.05 

1.0(f) 

2.13 ±0.19 

0.0(f); 6.48 ±0.11; 1.66 ±0.50 



abundance set to the solar value); //) an unabsorbed PL with 
photon index F ~ 2.1 ± 0.2; and //;) a narrow (cte ^ 0) gaus- 
sian Une ati? = 6.47±0.12 keV, with equivalent width EW -- 
1.7 keV. In Fig. 3 we show model D {top left, with unfolded in- 
strumental response) and the residuals of the fit {bottom left), as 



Fig. 3. The model spectrum C, unfolded of instrumental re- 
sponse {top left); and the fit residuals {bottom left). The 68%, 
90% and 99% probability contours for two interesting parame- 
ters for: (top right) the temperature, kT (in keV), and the chem- 
ical abundance, Z (in solar units), of the thermal component; 
(middle right) the intrinsic photon index, F, and the in situ HI 
column density, iVn, for the PL component; and (bottom right) 
for Fe-K line energy vs. normalization. 



well as the confidence countours of some physically meaning- 
ful parameter pairs (right). [If we allow for intrinsic absorption, 
none is deduced in the soft thermal component, while some is 
found (at > 90% confidence level) in the PL component, corre- 
sponding to TVh = (1.0 ± 0.8) X 10^^ cm^^. Intrinsic absorp- 
tion would cause the PL slope to steepen to F = 2.7 ± 0.5; see 
Fig.3-middle right.] 

The spatially integrated 0.5-10 keV spectra of local SBGs 
are well fitted by a model including a soft (kT < 1 keV) ther- 
mal emission, and a higher energy component consisting of ei- 
ther a ?J; 5 keV thermal emission, or a F '-^ 1.5 — 2 PL (e.g., 
Dahlem et al. 1998). By its continuum properties, therefore, the 
BeppoSAX spectrum of NGC 4666 seems to be a typical SBG 
spectrum. 



5. The Fe-Ka emission iine 
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Table 5. Spectral analysis of XMM-Newton data 



What makes the integrated spectrum of NGC 4666 quite out- 
standing among known SBG spectra is the huge (EW ~ 2 keV) 
emission feature clearly required by the BeppoSAX data at ~6.5 
keV (see Fig. 3, bottom right). In fact, in SBG spectra emis- 
sion between 6 and 7 keV is uncommon, inconspicuous (EW 
^ 0.3 keV), and centered at 6.7 keV (and hence it is identified 
with Ka emission from high-ionization iron - e.g., Persic et al. 
1998; Griffiths et al. 2000; Pietsch et al. 2001). 

As this emission is unresolved in the BeppoSAX data, we 
search for it in the archival XMM-Newton data. The superior 
spatial resolution of the latter allows us to identify the line as 
originating from the nuclear region, and clearly distinguish it 
from the extended continuum emission (see Fig.4). 

The nuclear location of the line, revealed by XMM-Newton, 
could be the clue to the presence of a strongly absorbed AGN 
in NGC 4666. If so, according to the Compton-thick model 
predictions we expect to observe a prominent (EW ^ 1 keV) 
line centered at 6.4 keV and superposed on a flat continuum 
underlying the line itself and cospatial with it. 

Spectral analysis of XMM-Newton data for the nuclear and 
circum-nuclear region does indeed show such a line {E ~ 
6.42 ± 0.02 keV with EW - 0.7 ± 0.2 keV), superimposed 
on a flat (F ^ 1.3) continuum (see Fig. 5 and Table 5). We 
identify the line as due to the Ka transition at 6.4 keV from 
low-ionization iron. 

To check for consistency of the Compton-thick hypoth- 
esis, we try a cold continuum reflection model (PEXRAV 
in the standard XSPEC spectral fitting package). This model 
represents the emission that would be expected from, e.g., 
a circum-nuclear torus that totally blocks the direct view of 
the primary AGN spectrum (i.e., the type-2 configuration in 
the unified model of AGNs). As commonly done for standard 
primary spectra of Syl and Sy2 galaxies (e.g. Turner et al. 
2000), we assume a PL with photon index F = 1.9 for both 
the primary (completely absorbed) component and the warm- 
scattered component (which is visible in the ~ 1-5 keV band). 
We also assume solar chemical abundances, full ionization of 
H and He, and a viewing angle of 80° (i.e., we take the torus 
axis to coincide with the spin axis). The resulting fit of this 
(astrophysically motivated) model is satisfactory (see Table 5), 
similar to the previous case of the (heuristic) single-PL model. 
The 2-10 keV flux of the best-fit cold-reflection component is 
-'^refl''^ ~ 4.8 X 10^^^ erg cm^^s^^, yielding an observed lumi- 
nosity of L^3^ ^ 4.0 X 10^^ erg s^^ (see Fig. 5). If we assume 
- as did Iwasawa et al. (1997) - that the intrinsic luminosity of 
the hidden AGN is a factor ^ 50 larger than the observed value, 
then L^^^> 2 x lO^i erg s^^. 

It should be emphasized that the extraction radius (^15") 
corresponds to a linear size of 1.9 kpc. Thus, this large circum- 
nuclear region includes a substantial fraction of the SB emis- 
sion. This may well explain why the continuum underlying 
the 6.4 keV line is somewhat steeper than predicted by a pure 
Compton-thick model (see Fig. 2). 



Model 



Parameters 



A: mekal+pl+gauss 
kT (keV) 

Z{Zq) 

r 

Line: cte; E(keV); EW (keV) 



x7/^ = 220.05/163 

0.33 ± 0.02 

1.0(f) 

1.31 ±0.05 

0.0 (f); 6.42 ± 0.05; 0.75 ± 0.25 



B: mekal+pl+gauss+refl. 

kT (keV) 

Z{Zq) 

F 

R (= refl./dir.) 

Line: o^\ E(keV); EW (keV) 



xV^^ = 236.35/163 

0.32 ±0.02 

1.0(f) 

1.9(f) 

6±2 

0.0 (f); 6.42 ± 0.05; 0.75 ± 0.25 



Fig.4. PN raw image (Ipix = 4 arcsec) in the 5-7 keV energy 
band {top), and in the 6.3-6.5 keV band {bottom). 



6. Origin of thie continuum 

6.1. Thermal component 

The thermal plasma component identified in the spatially in- 
tegrated BeppoSAX data agrees (within errors) with the ear- 
lier finding (based on ROSAT data) by Dahlem et al. (1998). 
Our BeppoSAX data are not sensitive enough to determine the 
chemical abundance, whereas the temperature is better con- 
strained to kT r^ 0.3 ±0.1 keV for a broad range of the abun- 
dance (see Fig. 3, upper right); also, the temperature is decou- 
pled from the slope of the PL component. 

Thermal emission is usually interpreted as arising from the 
interaction between the outgoing galactic wind and the ambi- 
ent ISM. The wind is hot but very tenuous, while the ISM is 
dense but cold, so none of these is likely to emit appreciably in 
the X-ray region. The boundary region between the two phases, 
however, does achieve a combination of temperature and den- 
sity suitable for sub-keV emission from a few percent of the 
wind mass (Strickland & Stevens 2000; Suchkov et al. 1996). 



Fig.5. Left: Spectrum from the nuclear(15"radius)region (see 
Fig.l, bottom). A very strong neutral Fe-Ka line at 6.40 keV, 
and the underlying cospatial flat continuum (approximated by 
a r ~ 1.3 PL model), are clearly visible. Right: 68%, 90% 
and 99% probability contours for the Fe-Ka line's energy ver- 
sus normalization. The emission line is clearly significant and 
spectrally localized. 



Single-temperature models {kT ^ 0.8 keV) and sometimes 
two- temperature models (with kT ~ 0.8 and ~ 0.3 keV) have 
been deduced from ASCA + ROSAT(e.g., Dahlem et al. 1998) 
and XMM-Newton (Franceschini et al. 2003) data of actively 
star-forming galaxies. ^ 

Analysis of archival XMM-Newton data shows that the soft 
(0.3-2 keV) emission is diffuse and essentially co-extended 
with the higher energy 2-10 keV emission (see Fig.l). Because 
the soft emission is largely dominated by the thermal compo- 
nent, which arises from the galactic winds that are driven by 
current SF activity, its spatial distribution suggests that the SB 
occurs throughout the disk in NGC 4666. 

6.2. PL component 

Persic & Rephaeli (2002; hereafter, PR02) have quantitatively 
assessed the roles of the various X-ray emission mechanisms 
in star-forming galaxies. They have used an equilibrium stellar- 
population synthesis model of the Galactic population of high- 
mass X-ray binaries (HMXBs) and low-mass X-ray binaries 
(LMXBs). The abundance of SN remnants (SNRs), from both 
Type II and Type la events, was also consistently estimated. 
From the literature PR02 derived typical spectra for these 
source classes. The spectral properties and relative abundances 
of the various classes of stellar sources allowed PR02 to calcu- 
late the composite X-ray spectrum arising from a stellar popu- 
lation of Galactic composition. 

PR02 determined that the ^2-15 keV emission is domi- 
nated by LMXBs; this is not surprising, given that the synthetic 
spectrum was calibrated by that of a normal, quietly star form- 
ing galaxy like ours. Due to the long times between the forma- 
tion and the onset of the X-ray phase (the optical companion 
is a subsolar main-sequence star) in LMXBs, their emission 
traces the higher, past Galactic SFR, and not the lower, current 



^ Part of the measured thermal emission may also come from OB 
stars, which are expected to be abundant in SB regions and whose 



SFR. Thus in order to obtain an estimate of the current SFR, 
the spectral contribution of LMXBs should not be included. 

Recent observations of nearby galaxies - mostly with 
Chandra - have revealed two new populations of X-ray point 
sources (XPs), not present in our Galaxy, with 2-10 keV lumi- 
nosities in the range 2 x 10^^ erg s~^ < L2-10 < 5 x 10"^^ 
erg s^^ (Very Luminous Sources: VLXs), and with L > 10"^^ 
erg s~^ (Ultra-Luminous Sources: ULXs) (Fabbiano & White 
2003) ^. In the 2-10 keV band relevant here, mean VLX and 
ULX spectra are best-fit by PLs with photon indexes of ~2 
and ~1.2, reminiscent of, respectively. Galactic Black-Hole- 
X-ray binaries (BHXBs) ^ (e.g., Zezas et al. 2002; Terashima 
& Wilson 2003; Foschini et al. 2002; see also Fabbiano & 
White 2003). VLXs and ULXs may describe different lumi- 
nosity/spectral states of BHXBs as a function of the mass infall 
rate onto the BH ^. The discovery of these sources has aug- 
mented the variety of known stellar endproducts that emit X- 
rays, and the PR02 method can easily incorporate such new 
classes (Persic & Rephaeli 2004). 

The population of VLX and ULX sources in a galaxy is 
more abundant the more active is star formation in that galaxy. 
Specifically, if the integral XP luminosity function (XPLF) is 
described as N{> L) ex L~°' (which corresponds to a differ- 
ential XPLF of the type ^ oc L"(i+")), based on a wealth 
of Chandra oservations it has been found that XPLFs are less 
steep in SFGs than in normal spirals and ellipticals. (For exam- 
ple, the two most nearby SB galaxies, have a ^ 0.5 (M 82) 
and a ^ 0.8 (NGC 253), while normal spirals have a ^ 1.2, 
and ellipticals have approximately a > 1.4; see Kilgard et al. 
2002.) So the integrated point source luminosity of SFGs is 
dominated by the highest-L sources (Colbert et al. 2003). If 
the XPLF in NGC 4666 is similar to those in other star-forming 
galaxies (e.g., a ^ 0.5 — 0.8; see Kilgard et al. 2002), then its 
spatially integrated 2-10 keV luminosity and spectrum will be 
dominated by VLXs and ULXs (unless there is a substantial 
diffuse hard component). The resulting spectral index, F, will 



^ The two quoted limits correspond to Eddington luminosities for, 
respectively, a 1.5 Mq neutron star and for a 8 Mq black hole (BH) 
which is the limiting BH mass obtainable via ordinary stellar evolu- 
tion. ULXs are called 'super-Eddington sources'. 

** In a few cases the optical counterparts of ULX have been identi- 
fied as O stars (Liu et al. 2002; Roberts et al. 2002b), thus suggesting 
a relation between ULXs and HMXBs, which of course are young 
objects that trace the current SFR. 

' Qualitatively, the inverse correlation between L2-10 and F (i.e., 
flatter spectra for higher luminosities) can be understood within the 
standard accretion model devised to explain the high Lx of interacting 
binaries. The increase in Lx is driven by an increase of M, which in 
turn means a piling up of material around the emission region that 
leads to a higher Compton scattering optical depth, and hence to a 
flatter spectrum. E.g., according to Sunyaev & Titarchuk (1980) and 
Shapiro et al. (1976), if a source of photons is embedded in a plasma 
cloud of optical depth r and temperature T, the escaping radiation has 
a high-energy photon index 



thermal emission is characterized by kT 
s"^ (seePR02). 



0.2 keV, L ~ W" erg 
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Fig. 6. 2-10 keV contours obtained with the MOSl overplotted 
on the optical image from the DSS (for details see Fig.l). The 
seven "sources" (i.e., maxima of emission) are identified by 
15"-radius circles. An eighth circle denotes an empty field that 
was used for background estimation. 



be 1.2^ f ^ 2, with the exact value depending on the XPLF 
index a and on F = r(L). Data with angular resolution at the 
arcsec level and sensitivity at the ^^10^^'' erg cm^^s^^ level 
(e.g., Chandra data) will resolve XPs in NGC 4666 down to 
-^2-10 ~ 5 X 10"^^ erg cm^^ s~^. The fluxes and - where pos- 
sible - the spectra obtained for these XPs would allow one to 
obtain some firmer clues on their nature. A full analysis of the 
maxima of emission seen in the 2-10 keV XMM-Newton map 
is beyond the scope of this paper. A preliminary analysis is, 
however, in order. 

From the 2-10 keV isophotes overlaid on the optical image 
of NGC 4666 (see Fig. 6), 7 emission maxima can be identified 
by visual inspection. In order to estimate the fluxes from these 
locations, we have constructed 7 circles, with radii of 15", cen- 
tered on the selected features (as to include >70% of the en- 
circled energy). An eighth, off-galaxy circle was used to esti- 
mate the background. Assuming for each "source" a flat spec- 
tral shape (F = 1.3) in order to obtain a (likely) upper limit to 
its 2-10 keV flux and a background as given by the off-galaxy 
empty field, the most obvious conclusion was that "sources" 
#1-4 qualified as ULXs and "sources" #5-7 qualified as VLXs. 
However, of such sources in Fig. 6, three quite clearly show 
structure: #1, #3 and #6. These emissions can be interpreted as 
resulting from (at least) 2, 3, and 2 sources each, respectively. 
In fact, concerning emission #1, we know that about half of 
the emission is due to "cold reflection" from a deeply absorbed 
type-2 AGN (L2-10 



0.40 X 10^" erg s ^, see section 



5): the remaining flux, corresponding to L2- 
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3.8 X 10^9 
erg cm^^ s^^, can be attributed to an ULX. As for emission #3, 
thatmaybecomposedof 3 VLXs withL2_io ~ lO'^^ergcm^^ 
s"'^ each. Finally, the elongated emission #6 may come from 2 
VLXs with £2-10 ^ 5 X 10'^* erg cm^^ s^^ each. Furthermore, 
emission #2 looks quite confused and its isophotes, which tend 
to become more triangular for deacreasing levels of emission, 
may suggest the presence of three unresolved sources: given 
the global level of emission measured for this feature, such 



/(1039 erg s-i)< 2). 



three sources could be (e.g.) 2 VLXs with ^2-10^ 10'^^ erg 
s^i each, and 1 ULX with L2-io?5 10^^ erg s"^ 

Based on these considerations, we tentatively identify the 
XPs hinted at by the 2-10 keV XMM-Newton image as being 
3 ULXs (with luminosities 2< L2-io/(10^^ erg s"i)< 4) and 

9 VLXs (with luminosities 0.5^ -^2-10 
Collectively, they contribute 1.6 x 10'*" erg cm^^ s 
2-10 keV band. Adding the estimated cold reflection from 
the hidden AGN, the total XP luminosity of NGC 4666 be- 
comes L^^io = 2 X 10"*" erg s"!, i.e. 80% of the total 2- 

10 keV luminosity observed with BeppoSAX (the latter being 
L2-W = 2.5 X 10^° erg cm^^ s^^, see Table 1). The remaining 
20%, missed by the XP luminosity computed from the XMM- 
Newton map as described above, might be genuinely diffuse 
(see below), or it could be due to unresolved sources of lower 
luminosity (i.e., with £2-10^ 5 x 10'^® erg s^*), or a combina- 
tion of the two. 

We can gain an insight into this issue by the following ar- 
gument. If the differential XPLF is of the type ^ ex L~(*+"\ 
our tentative identification of XPs suggests a ~ 0.5 (similar 
to what is observed in M 82, see Kilgard et al. 2002): such 
XPLF implies that the collective observed luminosity due to 
our tentative XPs amounts to 76% of the XP luminosity inte- 
grated down to 10"^* erg s~^ (i.e. down to the lower luminos- 
ity limit defining VLXs). This simple model implies that the 
XP luminosity, integrated between 10"^® erg s^^ and 4 x 10'^^ 
erg s-\ is Lf^io = (1/0.76 x 1.6 + 0.4 = 2.5) x 10''° erg 
s^^. This matches the integrated 2- lOkeV luminosity seen with 
BeppoSAX. What about the spectral shape? If we assume VLXs 
to have F = 2.0 - 2.5 and ULXs to have F = 1.3, our toy 
model implies F ~ 1.7 — 2.0, in reasonable agreement with the 
BeppoSAX result F = 2.1 ± 0.2. Therefore a toy model for the 
XP population, that involves ULXs and VLXs distributed ac- 
cording to an a = 0.5 XPLF plus a Compton-thick AGN, can 
bring BeppoSAX data and XMM-Newton data into reasonable 
agreement as far as both the luminosity and spectral properties 
are concerned. 

That the 2-10 keV emission of NGC 4666, once the 
AGN emission has been subtracted, is indeed mostly related 
to current SF activity can be checked as follows. From a 
sample of SB-powered ultra-luminous IR galaxies (ULIRGs; 
Franceschini et al. 2003) it can be deduced that the 2-10 keV 
to FIR emission ratio of pure SBs is 9 = log(L2-io/iFiR,) ^ 
-4.3 ± 0.3 (Persic et al. 2004). Upon subtraction of the AGN- 
related emission, NGC 4666 has ^ -4.0: so its disk 2-10 
keV emission is roughly at the level implied by its current SFR. 
This deduction seems reasonable because NGC 4666 exhibits 
the "smoking gun" proof of being a global SB: a large-scale su- 
pergalactic wind and a substantial radio halo. It should also be 
recalled that NGC 4666 is member of a small interacting group 
of galaxies (Garcia 1993); such a dense environment may well 
have tidally triggered a global SF activity in NGC 4666. 



6.2.1 . Compton scattered radiation 

As is well known, Compton scattering of (the radio- 
synchrotron emitting) relativistic electrons by the FIR and cos- 



mic microwave background (CMB) radiation fields results in a 
PL component with roughly the same photon index, F ^^ 1.8, 
as that of the radio emission (Schaaf et al. 1989; Rephaeli et al. 
1991; Dahlem et al. 1997; Sukumar et al. 1988). This closely 
resembles the observed 2-10 keV spectral profile. Compton 
emission may contribute ^ 20% of the observed 2-10 keV flux, 
as suggested by the following argument. [Formally, this is the 
level of 2-10 keV emission needed to supplement if?io (as 
measured by XMM-Newton) in order to reproduce L2-10 (as 
measured by BeppoSAX).] If the relativistic electrons and seed 
FIR photons are cospatial with a uniform magnetic field B, the 
Compton luminosity can be estimated from Lc = LsUph/UB, 
where Lg is the total synchrotron luminosity, C/ph is the energy 
density in the seed photon field, and Ub is the energy density in 
the magnetic field. Integrating the radio spectrum over the 0.01 - 
200 GHz range (the upper limit of the integration contributes 
very little to the integral), we get Lg ~ 5.1 x 10"^^ erg s^^. For 
a spherical region of radius r, the seed photon energy density 
is Uph — 3-Lboi/(47rr^c). Assuming L^oi — -^fir and r ~ 3 
kpc (see Dahlem et al. 1997), we have Uph ^ 2 x 10~^^ erg 
cm""^. If i? ~ 15 iJ,G in the disk where the FIR emission origi- 
nates (Dahlem et al. et al. 1997), then Ub = B^/{8tt) - 10"" 
erg cm~^. The Compton luminosity of NGC 4666, integrated 
over the entire range of upscattered energies is expected to be 
Lc ~ lO'"^ erg s^^, suggesting that FIR-Compton emission 
may account for a non-negligible fraction of the observed 2- 
10 keV emission, ^2-10 — 2.5 x 10^" erg s^^. However, this 
estimate depends sensitively on the mean value of B which is 
not known very well. The corresponding estimate of emission 
from Compton scattering off the CMB (mostly) in the halo is 
even less secure, due to substantial uncertainty in the mean halo 
value of B (even ignoring pertinent spatial considerations), but 
is likely to be weaker than that of the disk. 



7. Conclusions 

The X-ray spectrum of the nearly edge-on superwind galaxy 
NGC 4666 was measured in the 0.1-10 keV band with the 
LECS, MFCS, and PDS instruments aboard BeppoSAX, as well 
as in the 0.3-10 keV band with the EPIC instrument aboard 
XMM-Newton. As seen by BeppoSAX, the emission is not spa- 
tially resolved, and significant detection occurs only at energies 
< 10 keV. Emission measured by XMM-Newton is spatially re- 
solved as coming from the whole disk. 

SB and AGN activities coexist in NGC 4666. The former, 
which extends over most of the disk, is revealed, in the X-ray 
region, by (i) diffuse thermal emission associated with SN- 
powered galactic wind, and by (ii) high-luminosity XPs that 
are well-known indicators of enhanced star-formation activ- 
ity. It is also indicated by (Hi) an XPLF index a ^ 0.5, sug- 
gested here by a preliminary XP analysis of the XMM-Newton 
data and similar to those derived for other actively star-forming 
galaxies. The AGN activity is revealed by a prominent Ka line 
from "cold" iron at 6.40 keV and its relatively flat underlying 
continuum, both coming from the nuclear region, that originate 
from the reflection of the primary continuum by the cold inner 
wall of the circumnuclear torus. 



More specifically, our results can be summarized as fol- 
lows. 

(a) At energies ^ 1 keV the spectrum of NGC 4666 is dom- 
inated by a ^ 0.25 keV diffuse thermal plasma component 
(of unconstrained chemical abundance), distributed throughout 
the disk, which traces the interaction between the hot, tenuous 
wind elements and the cool, denser ISM in the disk. The distri- 
bution of this sub-keV plasma shows that the current SB phase 
involves the whole galaxy. 

(b) At energies ~^2-10 keV the integrated spectrum is dom- 
inated by a steep (F ~ 2) PL component. This emission is 
probably due to unresolved XPs distributed according to a lu- 
minosity function with (differential) index a ^ 0.5, and hav- 
ing L2-10 > 10'^^ erg s^^. These sources, which are likely to 
be accreting binaries with <8M0 black hole components, are 
known to dominate the XP luminosity of nearby star-forming 
galaxies and in the relevant energy range their spectra have 
F - 2 PL (for 0.1 < i2-io/(10^^ erg s^^) < 2) and F - 1.3 
PL (for L2-10 > 2 X 10^9 erg s'^). A F - 1.8 PL contribu- 
tion from Compton scattering of relativistic electrons (whose 
observed radio index is a^ ^ 0.8) by the ambient FIR photons 
may add a truly diffuse 2-10 keV emission. 

(c) The prominent (EW ^ 2 keV) emission line, seen but not re- 
solved in the BeppoSAX data, is identified in the XMM-Newton 
data as fluorescent neutral Fe-Ka at 6.40 keV originating from 
the nucleus of NGC 4666. This, together with the presence of 
a flat (F ^ 1.3) continuum in the same nuclear region that 
emits the line, suggests the existence of a strongly absorbed 
(Compton-thick) AGN, whose intrinsic 2-10 keV luminosity is 
L2_io^2xl04iergs-i. 

Compared with other composite AGN/SB galaxies, 
NGC 4666 shows an important difference. For example, in 
NGC 6240 and Arp 299 the AGNs have high intrinsic X-ray 
luminosities (^ 10'''^ — 10^'' erg s^^) and are seen according 
to the type-1 view (Vignati et al. 1999; Delia Ceca et al. 2002). 
Therefore in these galaxies the AGN dominates the emission 
at E^ 2 keV: the highly-absorbed primary continuum, trans- 
mitted through an obscuring cold medium, emerges from the 
absorber and dominates the emission at E^ 10 keV; and the 
unabsorbed primary continuum, scattered into the line-of-sight 
by a warm gas located outside the absorbing medium, dom- 
inates at 2keV^ E^ lOkeV. Only nt E ^ I keV can the 
SB be recognized from the thermal signatures of the wind 
components. (A similar situation is encountered in compos- 
ite AGN/SB-powered ULIRGs: Franceschini et al. 2003 and 
Persic et al. 2004). The situation in NGC 4666 is very differ- 
ent. Here the AGN has a much lower (estimated) intrinsic lumi- 
nosity (^2x 10''^ erg s^^) and is seen in type-2 orientation (the 
galaxy is edge-on). The combination of these two factors makes 
the AGN contribution to the 2-10 keV emission of NGC 4666 
quite small, ^15%, so the SB - far from being outshone by the 
AGN - dominates the integrated 0.3-10 keV continuum. 

The case that a low-luminosity type-2 AGN is identi- 
fied in the middle of the overwhelming glare of a SB makes 
NGC 4666 an observationally interesting and quite rare case. 
This occurrence owes to the simultaneous use of satellite 
data having widely different characteristics. The low-resolution 
BeppoSAX data have revealed the dominant continuum com- 
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ponent integrated over the disk, therefore giving information 
on the disk-wide SB; they have also revealed the existence 
of a prominent 6.4 keV line, apparently unrelated to the ob- 
served continuum. The higher-resolution XMM-Newton data 
have shown that the 6.40 keV ¥t-Ka line originates from 
the central region and, upon zooming onto such region, they 
have revealed a locally flat continuum, hence uncovering the 
existence of a Compton-thick type-2 AGN; they have further 
revealed a number of possible X-ray point sources that, in- 
tepreted as accretion-powered binary systems associated with 
active ongoing star formation, can naturally explain the inte- 
grated continuum observed by BeppoSAX. 
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